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This  paper  presents  an  attempt  to  estimate  the  parameters  of 
i  recursive  system.  The  technique  to  be  proposed  is  an  iterative 
■irnul at ion -estimation  procedure.  Thus  the  procedure  accounts  for 
he  fact  that  the  parameters  of  a  dynamic  system  are  being  es- 
imated . 

The  model  revolves  around  a  matrix  called  the  kedwood  Matrix 
:edwood[].  This  is  a  form  of  a  state  transition  matrix  whose 
tructure  is 
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where  m.  is  the  ratio  of  net  migration  of  i   year   olds 
to  the  total  popultaion  of  i  year  olds 

s.  is  the  probability  of  surviving  from  age   i   to   age 
i  +  1 

b.  is  the  number  of  expected  births  tor  eacn  member   or 
age  group  i 
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tfe  will  denote  the  matrix  in  (1)  by  P. 

The  recursive  population  projection  equation  is 

kt+i  -  pkt  <2> 

where  k   is  a  nxl  vector  whose  elements  represent  the   population 
at  time  t  of  the  n  different  age  groups. 

As  we  are  all  aware,  nowever,  both  social  and  economic  fac- 
tors can,  and  do,  affect  the  population.  Medical  factors  affect 
survival  rates.  Local  economic  conditions  may  warrant  a  change 
in  migration  habits.  Changing  per  capita  income  cnanges  the  op- 
portunity cost  of  having  a  child.  A  full  discussion  of  social 
and  economic  variable  affecting  population  may  be  found  in  Red- 
wood [  ]  . 

To  reflect  these  dependencies,  Redwood  regresses  the  birth 
rates  onto  various  sociological  and  economic  indicators.  He  did 
not  foresee  any  need  to  regress  the  survival  and  migration 
statistics  onto  the  time-varying  social  ana  economic  statistics. 
He  assumes  that  survival  rates  will  increase  steadily  over  time 
and  that  migration  rates  are  do  not  change  over  time.  I  will  as- 
sume constant  survival  rates. 
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Redwood  specifies  the  birth  rate  regressions: 

b*-4.i   =  ai   +  ai  At-  +  ai  ut-  +  ai  Lh  +  ai  h*  +  ai  W*-  +  ai  Dt- 
r.  +  ii     i0     i1  t     i2  c     i3  c     i4  u     i5  t     i6  ui 

(3) 


where 

A.  is  the  level  of  affluence 

U   is  the  unemployment  rate  at  time  t 

L.  is  the  female  labor  participation  rate  at  time  t 

W   is  a  dummy  variable,  1  =  war  year,  0  =  otherwise 

M   is  the  median  age  of  females  at  first  marriage. 

The  P  matrix  is  now  composed  of  time   varying   birth   parameters 

Denote  P  at  time  t  by  p    The  model  then  becomes 


kt+i  ■  ptkt  •  (4) 
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Estimation  Tecnnique 


It  is  not  the  interest  of  this  paper  derive  with  unbiased 
stimators  of  the  t^'s,  s_'s,  or  m^'s.  It  is  the  interest  of 
his  paper  compute  with  estimates  of  the  b.«s#  s-'s,  hl's,    such 

1        J        K 


hat 


T 

I    (k  -k  ) ' (k  -k  )  (5) 

t=l 


s  minimized,  wnere  k   is  the  nxl  observation  vector   of   popula 
ion,  and  kfc  is  evolved  by  the  difference  equation 

Kt+i  =  PKt  <6> 


vhere 


ke  =  ke  • 


The  reader  should  note  the  difference  between  (5)  ana  the  quanti 
ty  minimized  by  ordinary  least  squares  regression  on  lagged  vari 
ables,  which  would  be 


T 

* 


I    (k.-k. ) ' (k.-k. 
t  =  l      u       L 


where 


kt+i  ■  p  kt  • 
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Estimation  of  P  shall  be  done  through  a  method  known  as  sen- 
sitivity trajectory  analysis.   An  introduction  to  methods  of  sen- 
sitivity analysis  is  found  in  Perkins []  although  the  presentation 
:here  is  in  differential  equations. 

A  quick  survey  of  control  theory  tells  us  that  systems  can 
)e  mathematically  described  oy  state  equations.  State  equations 
lescribe  a  functional  relationship  between  the  current  endogenous 
ind  exogenous  variables  and  next  period's  endogenous  variables. 
jet  us  set  up  a  system  wnere  the  state  (endogenous)  variables  are 
issumed  to  be  of  the  form 

xk+i  ■  Axk  +  uuk  <7» 

/here  x  is  the  n  oraer  state,  (the  endogenous  variable),  u  is 
:he  m  order  exogenous  variable.  A  and  B  are  matrices  of  order 
ixn  and  nxm  respectively.  It  is  the  purpose  of  this  part  of  the 
>aper  to  identify  A  and  B;  that  is,  to  estinate  a  set  of  r  param- 
eters called  oC,  whose  value  bears  some  functional  relationship  on 
\  and  B,  say 

A  =  A  (ox) 

B  =  B(<X) 

)ur  state  equation  becomes 

xk  +  l  =  A^)xk  +  &(«)uk  .  (8) 
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,et 


X1  =  2_ 


oo^k 


(9) 


vhere  N^   is   a   nxl   vector   whicn   represents   the   sensitivity 


:  th 


:ra jectory  tor  the  change  in  the  parameter  <£.     (the  i    element  of 


*)  .  This  method  is  similar  to,  yet  easier  to  use  than  the  x.  's 
iiscussed  in  Box  and  Jenkins [, pp. 223] .  The  sensitivity  trajecto- 
ry approximates  the  peturbation  is  the  state  trajectory  due  to  a 
perturbation  in  a  parameter  of  the  state  equation. 


By  differentiating  (9)  with  respect  to  <£.  we  have 


ScTT  xk+l  "  ^k+1 

A(«  ]K   +  8ff7A<«   >Xk   +  5arD<   «   >uk 


(10) 


It  we  perturo  the  parameter  c^  Dy  the  amount  6c<1 ,   then   at   eacn 
time  period  k  the  perturbation  in  the  state  is  approximated  by 


xk  -  xk  +  6cc±xJ 


(11) 


Let  us  now  apply  the  above  to  tne  problem  at  hand.  This  will 
be  done  in  two  steps,  first  for  single  parameter  estimation,  then 
for  multiple  parameter  estimation. 

As  was  stated  in  (5) ,  we  desire  is  to  minimize  total  simula- 
tion error.   We  have 
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min      T 

*.      ^^t-tO'^ftt)  (12) 


</here    k      ana    k.    are   defined    as    in    (5). 
As    in    (11)    we   have 

kt    =    k"   +  6cCiX^    -  (13) 

*!?  and  ^J  are  the  state  and  the  corresponding  sensitivity  trajec- 
:ory  when  cC.  is  set  ot  the  nominal  value  or  ex" .  Trie  <Xl  is  some 
Jjfb-,  or  mk  parameter  in  (1).   Therelore  we  can  approximate  (12) 


Jx<kt  -  <k?  +zz*«0)'<*t  -  <k?  +  oij6^)) 


)r 


T  ok*.  ok, 

Jx<   et"  3g-6*i)'<et-  3S-6*i>  <14> 


vhere 


et  =  k     -  k£ 


The    first    order    condition    for    minimization   of    (14)    is 


T  Ok  ok 

0  =  "2t!1^(et-3^fei) 


Solving    for    6&i    we    have 
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GcC-  = 


—  t  e 

2  o^-  t 

t  =  l   i   

T  okt ,okt 


(15) 


We  perturb  cC  by  tnis  value,  tnat  is  we  redeLine  c<"  to  De 
!?  +  6cC^ .  The  the  procedure  is  repeated  for  the  new  value  of  the 
parameter.  The  state  and  sensitivity  trajectories  are  recomput- 
d.  Using  these  new  trajectories  another  perturbation  of  the 
>arameter  is  computed  by  (15).  This  procedure  is  continued  until 
:he  parameter  perturbations  become  sufficiently  small. 


-s 


As    in    (1C)    the  difference   equation    tor    tne    sensitivity    state 


d^"Nt  +  l 


i  ok. 

6<X         t  cxX. 


(lb 


vhere 


s^-P  is  a  matrix  of  zeroes  with  a  one  where  the  parameter  <£   ap- 


6k, 
5* 


=  0  . 


ears . 


The  proceoure  tor  the   simultaneous   estimation   of   several 
parameters  is  discussed  below. 


If  the  state  equation  is 


t  +  1  =  A<«  >kt 


Then  the  state  equation  tor  each   sensitivity   trajectory 
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.s 


'V  +  l  ~  5^"kt  +  l~A(3:)^t+5^-Akt 


(17) 


Given  initial  values  ior  each  element  ot  parameter  vector 
we  can  approximate  the  perturoation  ot  the  state  trajectory  by 
i  perturbation  of  the  parameter  vector  by 


1=1 


(lb) 


[12)     is    approximated    by 


nun 


1 


«      J^t-^'^t"^ 


(19) 


,et 


kt  -   k"t  "     2     fei*t 

1=1 


=   efc    -   l   6oCi>,j 


(20) 


=   efc    -    Lfc6c<      . 


(21) 


Then    (19)    becomes 


mm 


T 


c*     J    (et-Lt6oc>  V(et-Lt6oO 


(22) 


tfhere 
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Lt    =    I    ^1 »^2»^3» *  *  * '^sJ 


et  =  i*t-k?) 


6cc  =    [60^,  6c*2,  6c<3,  ...  ,6^]  ' 


vie   can    rewrite    (22)    as 


nun 


V  '  V 


(23) 


tfhere 


v    = 


e,    -  l16cc 
e3  -  L36c; 


eT  -  lt6oc 


Tne    rirst    oraer    condition    Lor    optimization    is 


0    =   X—v'v    = 
0    (X 


-2L'  (e    -    L6a:)     , 


(24 
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vhere 


L   = 


and 


e    = 


T 


4 


The    vectors    e   and    v   are    n'ixl,    while    L    is    a   nTxs    matrix.       (23)    be- 
comes 


L  '  e    =    L  '  L6cc 


(25) 


or 


6cx    =     (L'L)     1L' 


(26) 
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we  perturb  the  values  ot  the  parameters  by  the  amount  speci- 
fied by  (26).  With  these  new  values  01  the  parameters  we  recom- 
pute the  state  and  corresponding  sensitivity  trajectories.  We 
compute  the  new  perturbations  ot  tne  parameters  that  will  minim- 
ize (23)  by  (26).  This  procedure  is  continued  till  the  perturoa- 
tion  of  the  parameters  are  smaller  tnat  some  prespecir ied  value. 
By  (26)  we  see  that  the  parameter  perturbations  become  small  when 
the  error  becomes  orthogonal  to  the  sensitivity  trajectories. 
This  is  in  constrast  to  ordinary  least  squares  estimation  in 
which  the  parameters  are  cnosen  such  that  the  error  is  orthogonal 
to  the  predetermined  variables. 

If  we  were  to  estimate  according  to  (4),  instead  of  (2),  all 
that  changes  is 


a^-kt+i  ■  ssr^th  +  Ptsq  •  (27) 


In    tnis    case,    cC  does    not    include    the    b    *Sf    but    the   ai's    from    (3) 
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